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1 INTRODUCTION 



ABSTRACT 

Redshift space distortion (RSD) is a powerful way of measuring the growth of structure 
and testing General Relativity, but it is limited by cosmic variance and the degeneracy 
between galaxy bias h and the growth rate factor /. The cross-correlation of lensing 
shear with the galaxy density field can in principle measure 6 in a manner free from 
cosmic variance limits, breaking the f — b degeneracy and allowing inference of the 
matter power spectrum from the galaxy survey. We analyze the growth constraints 
from a realistic tomographic weak lensing photo-z survey combined with a spectro- 
scopic galaxy redshift survey over the same sky area. For sky coverage /sky = 0.5, 
analysis of the transverse modes measures b to 2-3% accuracy per Az = 0.1 bin at 
z < 1 when 10 galaxies arcmin"^ are measured in the lensing survey and all halos 
with M > Mmin = 10^^h~^ Mq have spectra. For the gravitational growth parame- 
ter parameter 7 (/ = ^m)^ combining the lensing information with RSD analysis of 
non-transverse modes yields accuracy (7(7) ~ 0.01. Adding lensing information to the 
RSD survey improves (7(7) by an amount equivalent to a 3x (10 x) increase in RSD 
survey area when the spectroscopic survey extends down to halo mass 10^^ '^ (lO^'') 
h~^MQ. We also find that the 17(7) of overlapping surveys is equivalent to that of 
surveys 1.5-2 x larger if they are separated on the sky. This gain is greatest when the 
spectroscopic mass threshold is 10^'^-10^''/i~^Mq, similar to LRG surveys. The gain of 
overlapping surveys is reduced for very deep or very shallow spectroscopic surveys, but 
any practical surveys are more powerful when overlapped than when separated. The 
gain of overlapped surveys is larger in the case when the primordial power spectrum 
normalization is uncertain by > 0.5%. 

gravi ty, for example, 7 is d ifferent from 7gr by more than 
20% (|Linder fc Cahnll2007l ). 



Measurement of the linear growth of structure of the Uni- 
verse is essential in that the growth history reflects the 
natu re of dark energy and the underlying gravity model 
(e.g. lYamamoto et al. l2010t ). i.e. whether dark energy is a 
cosmological constant, or is evolving with time, or if Gen- 
eral Relativity (GR) is the correct gravity model that gov- 
erns the evolution of the Universe. In the linear regime of 
GR, the growth of perturbations is scale independent. It 
can be parameterized as the linear growth function G, with 
P{z) = G^{z)PcMB, where P{z) and Pcmb are the mat- 
ter density power spectra at redshift z and at the epoch of 
recombination, respectively. G{z) carries information about 
the amount of dark energy and dark matter. The growth 
rate factor f = a'"'^ , with a being the scale factor, is an- 
other quantify of interest: / can be well approximated as 
/ = QZi, with 7 in a narrow range near 0.55, for a w ide vari- 
ety of dark-energy mod e ls in General Relativ ity (|Peeblesl 
1 19801 : iLahav et al] Il99ll : iLinder fc Cahiil |2007| ). A precise 
measure of 7 therefore enables one to distinguish GR from 
alternative gravity models. In a braneworld type of modified 



Precise measurements of G and / constrain dark en- 
ergy and gravity, and Redshift Space Distortion (RSD) has 
been shown to be a powerful approach to perform this mea- 
surement (e.g.LK aiscr'l987'; 'Colc ct al."l994'; 'Hamilto n et al 



200c : IPeacock ct al.. .2001.: .Scoccimarro, .2004: Guzzo et al. 
20081 : ICabre fc Gaztanagal HoogI : Ipiake et al.1 I2OIII ). This 



RSD measurement is, however, only precise in the linear 
regime. At late epochs of the Universe, the linear regime 
(of the velocity field in particular,) is confined to very large 
scales, k < O.l/i^^Mpc^^. On these scales, the measurement 
is usually limited by sample variance, or cosmic variance — 
we do not have many independent perturbation modes for 
the measurement because of the finite survey volume ob- 
servable in a given epoch. 

Using multiple tracers of the density field, one 
can in principle evade sampling variance, and mea- 
sure the linear growth of s t ructu r e with unbounded 
accuracy (iMcDonald fc SeliakI [20091: IWhite et al.1 l2009l : 
iGil-Marm et al.l I2OIOI : iBernstein fc Gail I2OIII ). The great 
benefit of multiple tracers is not realized, however, if only the 
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clustering of galaxies is measured because of the following; 
using RSD, one can only measure j3 — f/b and the prod- 
uct fG. Without any prior knowledge of the galaxy bias, 
one can not constrain / o r G in dependently. It has been 
shown by iBernstein fc Gail (|201ll )[BCll] that prior knowl- 
edge on galaxy bias significantly improves the constraint on 
the growth of structure in the case of single survey redshift 
bin. 

In principle, galaxy bias can be measured by cross- 
correlating weak gravitatio nal lensing convergence with 
galaxy clustering (|Penll20()3 ). This bias measurement is free 
of sample variance in the sense that the bias errors from a 
survey of a fixed number of modes can be reduced without 
limit as the lensing measurement noise and galaxy shot noise 
are decreeised. This enable us to use large-scale modes that 
are well in the linear regime for high-accuracy measurement 
(see ECU). Combining a weak lensing survey (using pho- 
tometric redshifts for source galaxies) with a spectroscopic 
redshift survey of lens galaxies can serve this purpose per- 
fectly: the bias of the spectroscopic galaxies is measured by 
cross-correlation with the lensing signal in transverse modes, 
while the RSD analysis of the spectroscopic sample is con- 
ducted using non-transverse modes over the same volume. 

Current and future large surveys are making possible 
the combination of spectroscopic redshift-space and lensing 
maps over common volumes. For example, the footprints of 
the upcoming Dark Energy Survey (DESj^may overlap with 
that of an extended Baryon Oscillation Spectroscopic Survey 
(BOSS^ survey near the equa tor. The future Euclid space 
telescope (|Laureiis et al.|[201lf ) is designed to take spectra 
and images of galaxies at the same time. 

In this work, we will explore the potential improvements 
in constraint of growth of structure and gravity from over- 
lapping RSD and lensing surveys. This is an extension of 
BCll, where we consider the case of one single RSD redshift 
bin with an arbitrarily assigned prior on galaxy bias. In this 
work, we will consider the more realistic case of tomography 
using spectroscopic (RSD) and photometric (lensing) sur- 
veys covering common sky area, with both types of tracers 
divided into as many as 20 redshift bins. We will explore how 
well galaxy bias can be measured using the cross-correlation 
of galaxy shear and galaxy clustering in this realistic joint 
tomographic survey. The basic scheme is: 

(i) Conduct a galaxy redshift survey and a weak lensing 
(photo-z) survey over the same volume of the Universe. Split 
both galaxy samples into redshift bins. 

(ii) Optimally weight galaxies in each bin of the redshift 
survey to pro duce a mass density estimator with m inimal 
stochasticity (|Hamaus et al.l I2OI0I : ICai et al.l boilf ). Mea- 
sure the 2-point shear-shear (from the lensing survey), and 
density-density (from the spectroscopic survey) correlations 
and the shear-density cross correlations between all z-bins. 
Using these measurements of the covariance in transverse 
modes, constrain the bias b of the spectroscopic galaxy den- 
sity estimator and the mass power spectrum Pm in each 
redshift bin. 

(iii) In each z-bin, split galaxies from the redshift survey 



into difi'er ent bias bins. Perform multiple-tracer RSD mea - 
surement (|McDonald fc SeliakllioogI : IBernstein fc Caill201ll ) 
using the redshift-space density field of the binned galaxies. 
The b and Pm constraints derived from transverse modes 
in step (ii) are incorporated to break the f — b degeneracy 



inherent to RSD, so that separate constraints on G and / 
can be achieved. Throughout the paper, we will suppress the 
latin index denoting redshift in equations that involve only 
a single redshift bin, such as the RSD Fisher matrix. We use 
Greek indices for galaxy bias bins. 

We use the Fisher matrix method to forecast growth 
constraints resulting from a model survey consisting of spec- 
troscopic and lensing surveys covering a common /sky = 0.5 
of the sky, reaching the depth of z = 2. We split both sam- 
ples into 20 z-bins of width Az = 0.1. We employ the halo 
model for our survey model, assuming that each halo above 
mass Afmin hosts one spectroscopic target galaxy. We set up 
our forecast methodology for lensing tomography in section 
[2] and for multi-tracer RSD in section (3] We summarize our 
numerical results in sectional and conclude in section [S] 

Unless noted otherwise, we assume a fiducial flat ACDM 
cosmology with the following parameters: Qm.=0-272, = 
0.728, Ofe = 0.0455, gg = 0.807, Us = 0.961, Ho = 70.2. 
(|Komatsu et al.ll201l] ') 



2 GALAXY-SHEAR CROSS-CORRELATION 
2.1 Weak lensing tomography 

Weak gravitational lensing of background galaxies is a pow- 
erful way to measure the projected mass density of the fore- 
ground. It is free from galaxy bias and can be used to mea- 
sure galaxy bias when cross-correlating with the galaxy den- 
sity fleld. Source galaxies split into different tomographic 
bins enable us to probe the mass density at different epochs 
of the Universe. 

For the ith 2-bin of source galaxies, the observable of 
weak lensing is the distortion of those galaxy images, or cos- 
mic shear, which is induced by foreground large-scale grav- 
itational potential. From the cosmic shear one can infer the 
convergence k, which is a weighted projection of the 3-D 
mass density of the foreground: 



2c2 



'"xw^.(x)^rfx (1) 

where x is the comoving radial distance, a is the scale factor 
of the Universe, 5 is the 3-D matter density contrast. The 
lensing weight function is 

dz Xs 



w^ix) = - 



ni{z) 



X 



dXs Xs 



dXs, 



(2) 



where fii is the number of galaxies in the ith redshift bin, 
distributed as ni{z). xn is the horizon distance. We assume 
a total source redshift distribution of the form 



n{z) = noz exp(— z/zo), 



(3) 



^ http:/ /www. darkcncrgysurvcy.org 
2 http://cosmology.lbl.gov/BOSS 



where zq — 0.45 is chosen to fit with the predicted Eu- 
clid survey's median redshift, and no is chosen such that 
/ n{z)dz = A'^ions, the total density of lens source galaxies 
per steradian. In section |4l we will examine results for a 
wide range of A'lcns- 
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The weak lensing signal is detectable only in statis- 
tics of large source-galaxy ensembles, e.g. via the two-point 
correlation function or its Fourier space counterpart, the 
shear pow e r spec trum, or higher order correlation functions. 
iBernsteinI (|2009l ') gives a framework for two-point analy- 
sis of weak lens ing survey data. We follow the notation of 
IBernsteinI (|2009l ') for the Fisher matrix from lensing tomog- 
raphy. We will work in the Fo urier domain , and employ 
the Limber approximation (e.g. lLimberlll954l : iKaiseij [19921 : 
IHuI|2000| : I Verde et al]|2000l : ICai et al.ll2009l ). assuming that 
there is no correlation between S in different redshift bins 
nor between different spherical harmonics. We also assume 
that within each redshift bin, ni{z) is a Dirac delta function 
at Zi, and that the projected mass fluctuations Si within bin 
i can be treated as a single lens deflection screen at Zi. Un- 
der these assumptions, for a given spherical harmonic, the 
convergence of the ith source galaxy bin is just the weighted 
sum of the mass density of all the redshift bins in front of 
the ith bin (i increases with redshift): 



(4) 



with (eiEj) — S('^a^/ni the variance of lensing shear noise. 
We take CTj = 0.22 throughout our calculation. Aik ~ 
, Di is the comoving angular diameter distance to 
redshift Zi, and Ft = |i?o^^m(l -I- Zk)- 



2.2 Covariance matrix for lensing and galaxy 
density 

When combining lensing tomography with a galaxy redshift 
survey over the same volume, we assume that the spectro- 
scopic galaxies will be split into z bins matching the source 
bins. A projected density estimator Sg will be produced in 
each bin using some weighted combination of the spectro- 
scopic galaxies. These projected density estimates are es- 
sentially the transverse modes of the RSD measurement in 
section (3] Each galaxy is given an optimal weight as de- 
scribed in section 12.41 The spectroscopic galaxies can have 
a different selection function from the lensing source galax- 
ies. In cases where spectroscopic galaxies are not available 
(such as when we consider non-overlapping RSD and lensing 
surveys), we will assume that the Sg measurement is made 
using galaxies with photometric redshift assignments from 
the lensing survey's imaging data. 

These measurements will be made for each mode trans- 
verse to the line of sight, indexed by spherical harmonic I: 

(i) Cij^ip) — the (cross-) power spectrum of the lensing 
convergence at (and between) different redshift slices from 
the lensing survey. 

(ii) Cfj{l) — the power spectra at each redshift slice of 
the projected galaxy density estimator formed from the 
weighted spectroscopic galaxy survey (or photo-2; sample). 
We assume no correlation between densities of distinct red- 
shift slices, following the Limber approximation and ignoring 
magnification biases and redshift mis-assignments, so C^^ is 
diagonal. 

(iii) Cf!^(l) — the shear-galaxy cross-spectra between dif- 
ferent redshift slices. Galaxy density will only correlate with 
shear in the background, so C®" is a triangular matrix. 



More specifically, the measurements can be expressed 

as: 

min{i,j } — 1 

CiJ^H) = AikAjkAxkPk{l)Fi +aUij/ni 

ci'ii) = A^'Axr'p.a/A)6?5,,+M(o 

C^;{1) = A,,D-'n{l/D,)hF U<j (5) 

where Ui is the number density of lens source galaxies per 
steradian at the ith redshift bin; Pi(k) = Gi PcMsik) is the 
3-D mass power spectrum at the ith redshift slice, with Gi 
being the linear growth function at Zi; I = k/ Di is the angu- 
lar wavenumber; and bi and Ni {I) are the scale-independent 
bias and stochastic noise power, respectively, of the weighted 
spectroscopic galaxy density estimator at Zi. Note that we 
ignore complications from intrinsic alignments of galaxies, 
photometric redshift errors, and other lensing measurement 
systematic errors. 

The fiducial value of the noise term Ni{l) is taken from 
a fiducial stochasticity Ei{l) of the galaxy density estima- 
tor: N!"'{1) = g»(Z)^A~'Axr'fi (^/A)b?<5», . Our model for 
the fiducial Ei{l) is taken from [Cai et al.l (|201ll ) and de- 
scribed in section [2.41 The fiducial value oi Mi{l) is taken 
from this model, but Ni{l) is still treated as an independent 
free parameter of C^f{l), and will be marginalized over. Note 
that most analyses (including our own RSD Fisher matrix) 
assume that the stochastic power is known a priori to be 
given by the Poisson formula. We find in this lensing anal- 
ysis that a strong prior knowledge of M can substantially 
infiuence the final growth constraints, so we adopt a weak 
conservative prior quantified below. 

For each spherical harmonic I, the full covariance matrix 
for the lensing and density measurements is 



C = 



which is a 40 X 40 matrix for our 20 redshift slices Zi 
0.1, 0.2... 2.0. 



2.3 Fisher matrix of the cross-correlation 

The free parameters of the model for the lensing Fisher ma- 
trix are: 

(i) The amplitude of the mass power spectrum at dif- 
ferent redshift Pi{l/Di), which is in turn a function of 
only the linear growth function Gi at redshift Zi, Piil) = 
GiPcMB{l/Di). In practice we use the parameter pi = 
In P,(0 = 2 In G, -I- const. 

(ii) The bias bi of the weighted spectroscopic galaxy den- 
sity. It is related to the biases bia of individual galaxy bins 
used in the RSD analysis by bi = X^a^ia^'d with the 
weights Wia as assigned in Section 12.41 How galaxies are 
made into different bias bins is detailed in Section [3] 

(iii) The noise Mi{l) in the galaxy-galaxy clustering mea- 
surement in transverse modes. 

We fix all parameters except for the 20 Pi, 20 bi and 20 A/i — 
note that this includes taking the cosmological distances Di 
as known. There are no lensing sources behind the highest 
redshift bin, so the last bias parameter is unconstrained. We 
therefore drop the rows and columns of the Fisher matrix 
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for the h, P, and A/" parameters of the highest redshift bin, 
leaving 57 x 57 elements at each multipole 



^'ElsW = ^Tr [c(0-'c(0,pC(/)-^c(0, 



(6) 



where p,q £ {Pi, bi,Mi}. The above equation holds true be- 
cause < K >= for the whole sky. For the first bin, there 
is no constraint from lensing, so the first row and column 
of are zero. We add priors on the Mi parameters to 
indicate uncertainties propotional to the fiducial values: 



iJ\fJ\f , prior 



(0 = S^: 



(7) 



Ni is the number of I bins that we use. This scaling produces 
a prior such that the mean A/i over all I bins in known to 
accuracy 2aMf'^. We choose the very weak prior q = 50 for 
our calculation. 

An example of the lensing Fisher matrix is shown in 
the left panel of Figure [T] for a single mode at I — 30. We 
find the Fisher matrix is close to block-diagonal, i.e. b's, P's 
and Af's at distinct z are only weakly correlated. Lowering 
the fiducial stochasticity A/i makes the Fisher matrix more 
diagonal. 



2.4 Sources of noise in the lensing measurement 

In the determination of the 6, using the cross-correlation 
of lensing and galaxy surveys, both the shear measurement 
noise and stochasticity between the tracer and the mass field 
serve as sources of error. We can rewrite Equation (|4]) as 



.(0 = ^A,fcAxfcFfc 



fc=i 



ekil) +£• 



- Si{l) + €i, 



(8) 
(9) 



with the galaxy overdensity 5l{l)=bk5k{l) + &k{l), where Ck 
is the stochastic component. The lensing observable hence 
has two indistinguishable stochastic components that are 
not properly traced by the galaxies — its shear measurement 
noise plus the total convergence from the mass fluctua- 
tions Si{l) = J2k^i ^ik^XkFkBkil)- Together these degrade 
the constraint on the mean bias, and also affect the con- 
straint on the mass power since b and P are strongly corre- 
lated (see the Fisher matrix of b and P at the middle panel of 
figure[T]). We will investigate in Section l47l1 how the choices of 
lensing source density Aliens and spectroscopic depth Mmin, 
which set these two noise levels, affect the constraint on b 
and P, and further affect constraints on the growth. 

Since the stochasticity of the galaxy density has been 
shown to be a limit for th e precision of weak lensing con- 
straints on the bias (e.g. iPenI |2004 ). we have incentive 
to reduce the stochasticity below the commonly assumed 
Poisson level. Sub-Poisson stochastici ty has been demon- 
strated in A-body simulation s , (e.g . iBonoli fc PenI l2009l : 
iHamaus et al.l |201G|: [Cai et all mm . Here we follow the 
method of ICai et al.l (|201ll ). hereafter CBS, for mmmiizme 
the stochasticity of a mass estimator from a weighted com- 
bination of halos. 

The optimal weight Wopt of each halo is a function 
of its mass and of the minimum mass J\/min of halo in- 
cluded in the survey. The resulting stochasticity between 



the weighted halo field and the mass field -Eopt- Explicit 
expressions from CBS for Wopt and iJopt are given in the 
Appendix. With this definition of -Eopt, the stochastic com- 
ponents of the galaxy density clustering can be written as 
< e{l)e'{l) >= _Bopt(0/[l - -Eopt(0]- This expression is a 
function of redshift, but for simplicity, we drop the latin 
index denoting redshift here. 

In this work, we use the CBS halo model description 
of -Eopt to produce the fiducial value of stochastic power 
Mi{l). In principle, both uiopt and -Eopt are functions of the 
Fourier wave number fc, the minimal halo mass of the cata- 
logue -Mmin and redshift z. However, since we find that lUopt 
depends very weakly on k in the linear regime, we will just 
adopt the Wopt for fc = 0.01/iMpc~^ at each Mmin and z. 
CBS shows how Eopt drops with Mmin', therefore, a deeper 
spectroscopic redshift survey targeting galaxies hosted by 
lower-mass halos leads to a higher-precision measure of the 
galaxy bias when cross-correlated with lensing. 

In using the CBS model for our fiducial value of E, we 
are assuming that the host halo mass of each spectroscopic 
target is known by some means, and that the spectroscopic 
targets are complete to the limiting halo mass. This assump- 
tion might be somewhat strong but it is not impossible for a 
real survey. One can imagine using the relatively deep lens- 
ing image survey to resolve satellite galaxies that are hosted 
by each spectroscopic galaxy's halo. The number of satellite 
galaxies could then be used to estimate the halo mass. We 
need a spectroscopic redshift of only the central galaxy of 
each halo. 

When we are forecasting scenarios in which there is no 
overlapping spectroscopic survey for lensing data at a given 
z, we assume that photometric redshift maps can produce 
a Sg density estimator with fiducial stochasticity Ei = 0.5. 
Note that the bias b for the photo-2 population can not be 
used in this case to constrain the biases of the spectroscopic 
population because of different selection functions. 



2.5 Summation over modes 

In the lensing Fisher matrix, we marginalize over all noise 
parameters A/i to have (19fe's -|- 19P's)^ left. We also 
marginalized over all those b's and P's of non-linear modes. 
In marginalizing over non-linear modes, we will retain more 
modes at high z, since for fixed I, the physical scale is larger 
at high z. Furthermore, non-linearity develops on smaller 
scales at higher redshift. Since accurate predictions of red- 
shift distortions will likely be available only in the linear 
regime, we do not use non-linear modes in our measure- 
ments. 

To separate linear and non-linear modes, we first as- 
sume that at z = 0.5, the linear modes have k < fc^ax ~ 
0.1/iMpc^^. We then compute the variance a'^{R!^f^,z — 
0.5) = ^ J k^P{k)W^{kR'^t)dk smoothed by a spheri- 
cal top-hat window function W{x) — 3[sin(2;) — a; cos(a;)]/a::'^ 



with the radius 
of -R^in at all 
a^iR'^Lz = 0.5) 



^min 

other 



We choose the scale 
so that a^{R^i^,z) — 
. = 27r/.RS 



— 27r/fcn^ax- 

redshifts 
We obtain fc^a: 
k^^^D{z) for each z. 

After marginalizing over non-linear b's and P's at each 
I, we sum the Fisher matrices for linear-regime parameters 
over all modes with Zmin < I < 'max, with Zmin = 10 for all 



loq[Fisher,„,(l = 30)] 

-7.0 -T.O -3. 



Lensing + redshift surveys 5 




Figure 1. Examples of Fisher matrices for 20 redshift bins for weak lensing source density Aliens = lOarcmin^ and spectrosopic survey 
depth M^in = 10^^ Mq; each survey covers f^]^y = 0.5. Left: Fisher matrix Flchs from joint lensing/galaxy measurements on transverse 
modes at I = 30, with parameters for bias 5, mass power P, and galaxy stochastic power Af at each z; Middle: Flchs for linear 5 and P 
after marginalizing over all b's and P's of non-linear modes, summing over all I, and marginalizing all A^'s; Right: Fisher matrix about / 
and P after summing Flchs + Frsd, after marginalizing over all biases. Each matrix block contains parameters from low z in lower-left 
to high z in upper right, as labeled, omitting the highest z bin which is unconstrained. Fisher matrices use logarithms of each parameter 
so that fractional errors are represented. The fiducial value of stochastic powers Af are estimated from halo mode described in section 
12.41 and a weak prior is applied. Note that correlations between distinct redshifts are always quite weak. 



redshift s: 

'max 

F£l,(/) = E /=^^(2' + (10) 

i = 10 

An example of the final lensing Fisher matrix is shown 
in the middle panel of Figure [1] While the constraints on 
bi and Pj are highly correlated for i = j , the corrrelations 
among b and P at distinct redshifts are very weak, and we 
can consider the experiment to give essentially independent 
results at every redshift bin. 



3 MULTI-TRACER REDSHIFT SPACE 
DISTORTION 

In this section, we review the basic idea of using redshift 
space distortion (RSD) to measure the growth of structure. 
This will be implemented in a spectroscopic survey. For each 
redshift shell, galaxies will be made into multiple bins of 
their parent halos' masses. Our measurements will include 
the redshift space power spectra of each sub-sample, and the 
covariance of all those galaxy bins. 

Each redshift shell of the spectroscopic survey will have 
an independent Fisher matrix. 

In the linear regime, galaxy overdensity 5" seen in red- 
shift space will be boosted relative to the matter overdensity 
S due to the large-scale inflow bulk motion of galaxies. The 
flrst-order large-scale peculiar velocity is related through 
the continuity equation to the linear growth rate factor 
/ = g . The redshift-space galaxy clustering therefore 
encodes information on the growth of structure. In Fourier 
space, iKaiscr. (1983) derives 

SUk) = ib^ + f^^)5{k) + e^, (11) 

with ec the stochastic portion of the galaxy density with 
(eaS) = 0. For the RSD analysis we assume a diagonal 



stochasticity matrix, (ea^a) = S^^/ric, i.e. noise in distinct 
galaxy bins is uncorrelated. be is the bias of the ath galaxy 
bin, and /i is the cosine of the angle between the k vector 
and the line of sight. 

Following BCll, the covariance of the multi-tracer RSD 
measurement is: 

C.,8(k) = Cov(5J(k),5^(k)) (12) 
= (b^ H- f^^^){bfi /Ai^)G^PcA/s(k) +£■^,3, 
£ap = (eQE/j)- (13) 

The free parameters in this measurement are: the biases 
be of the galaxy bins, assumed to be scale independent; 
the growth rate /, the linear mass power spectrum P — 
G^PcMB, where G is the linear growth function and Pcmb 
is the power spectrum at the epoch of recombination. For 
each mode at each re dshift bin, we have the Fisher matrix 
(|Tegmark et al.lll997l ) of RSD: 

FrsD(k) = ^Tr [C(k)-^C(k),,C(k)-^C(k),,] (14) 

where p, g £ {/, P, bi, b2...bMh}- We assume that galaxies are 
binned by the mass of their parent halos, use three log mass 
bins for each decade of mass (A'^;, = 15 mass bins for the 
case of Mmin = 10" Mq). The size of the Fisher matrix is 
{Nb + 2) X (Ai'i,-|-2). This multi-tracer RSD method improves 
the constraint of fG by a factor of up to 6.4 compared to 
the standard RSD method, where all galaxies are placed in 
one single bias bin (BCll). Without any prior knowledge 
of galaxy bias, neither method can constrain / or G alone, 
only the product fG. 

Nonlinear k modes will not be used in our analysis, 
applying the criteria from section [23] We sum over modes 
within fcmin < k < fcniax, whcrc fcmin = lininD{z), to yield a 
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total RSD Fisher matrix for our redshift bin: 
F^lD(k) = 7^ E A/, f; 47rfc2F^li3(fc)Afc. (15) 

Here V is the surveyed volume within the redshift bin under 
consideration. When integrating the RSD matrix over fj,, we 
are careful to remove a section around = representing 
the number of transverse modes used in constructing Flohs 
for the same bin of redshift and k — l/D. This avoids double- 
counting the information in the transverse modes of the 
spectroscopic survey if we are analyzing overlapping RSD 
and lensing surveys. For non-overlapping surveys, we do not 
exclude the fj, — modes from the RSD information. 

3.1 Combining RSD and lensing Fisher matrices 

Note that the RSD Fisher matrix is degenerate in the /- 
b direction. Combining RSD measurement with a lensing 
survey can break the degeneracy between b and /, and hence 
yield a tighter constraint on 7. There may, however, be non- 
zero covariance between box P values at different redshifts in 
the lensing measurement. So when combining the constraints 
from lensing with those from RSD, the parameters in each 
z bin can not be treated independently. We will need to 
create a large joint Fisher matrix for biases, /, and G at all 
redshifts, so we concatenate F^gj-, from all 19 redshifts into 
a single block-diagonal RSD matrix. 

We have to convert the lensing and RSD Fisher ma- 
trices to encompass a common set of parameters, then sum 
them, being careful not to double-count information. The fi- 
nal Fisher matrix will contain entries for the A'';, bias values 
bia, plus the growth and growth rate d and fi at each of 
the 19 measurable redshift bins, giving a final dimension of 
19 X (A^i, -|- 2) . The matrix is nearly block-diagonal with iso- 
lated redshift blocks, because Frsd is completely decoupled 
between redshift bins, and Flghs is nearly so. We retain the 
full matrix, however, for completeness. 

The lensing Fisher matrix elements for galaxy bias refer 
to the weighted mean bias hi = Wiabia for each redshift. 
We can convert the constraint on the weighted mean bias 
into a joint constraint on the individual bins' biases using 
the known weights Wia from the Appendix. 



4 RESULTS 

In this section, we compare the constraints on the growth of 
structure from having a lensing photo- 2; survey, a spectro- 
scopic redshift survey, and the combination of them. We will 
also investigate the case of having the two surveys over sep- 
arate volumes. We will explore how the results may depend 
on the depth of the spectroscopic survey and the photo-z 
survey. For the spec-z sample, we usually label the survey 
depth as the minimal halo mass Mmin , since we assume that 
the spectroscopic targets are the central galaxies of all ha- 
los with M > Mmin. Figure [2] plots the space density of 
targets vs Mmin at a few nominal redshifts, plus the to- 
tal projected sky density of targets vs Mmin- For example, 
Mmin ~ 10^^/i~^Mq corresponds to galaxies of Milky Way 
size or larger, with a space density of ~ 10~^'^ h^Mpc~^ 
and sky density ~ lOarcmin"'^; having Afmin ~ 10^^ Mq is 




Figure 2. Galaxy number density rispcc in the spectroscopic sur- 
vey versus tlie minimum halo mass Mmin a^t three different red- 
shifts estimated from halo model. We assume one spectroscopic 
target per halo. Black dashed line shows the projected galaxy 
number density versus Mmin at < z < 2. 

like a Luminous Red Galaxy (LRG) sample, with a space 
density of ~ 10~^'^ h^Mpc~^ and sky density ~ larcmin"'^; 
and Mmin = 10^*/i~^A/q is a rich cluster survey, with a 
space density of ~ 10~^/i'^Mpc~^ at z < 1 and sky density 
~ 0.01 arcmin"^. Keep in mind that the survey targeting 
M > Mmin generally yields the best possible cosmological 
constraints for a given target density. 

We will assume that the primordial CMB power spec- 
trum is known exactly, unless specified otherwise. We will 
show that for most cases, knowing Pcmb to 0.5% gives about 
the same growth constraints as fixing it. 

4.1 Lensing constraint on b and G 

We first examine the constraints on b and P at different 
z bins from FLcns, the joint analysis of lensing and galaxy 
density surveys in purely transverse modes. Assuming the 
primordial power spectrum Pcmb is known, measuring P is 
the same as measuring the linear growth function G. 

The dotted lines in Figure [3] plot the Fisher uncertain- 
ties in hi and Gi (the equivalent of y/P) vs redshift Zi. 
Each plotted point gives errors after marginalization over 
all other parameters. We find the measurements of galaxy 
mean bias and G reach percent-level accuracy over a large 
range of redshifts for Mmin = 10^^ h~^MQ (left panel), and 
sub-percent accuracy when galaxy stochasticity is lower with 
Afmin = W^'^ Mq (right panel). The constraint is better 
at low redshift, easily understood since higher-redshift lenses 
have fewer background galaxies to lens and hence higher ef- 
fective shape noise in the lensing measurement. 

The number of available linear modes increases rapidly 
at higher z, which should in principle cause constraints to 
improve with redshift. While we see this behavior at z < 0.5, 
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Figure 3. Lensing constraints on b and G (dotted lines), RSD constraint on fG (dashed lines) and Icnsing+RSD constraint on / and 
G (solid lines) in each redshift bin of width Az = 0.1 are plotted vs redshift, for the case of Anions = lOarcmin^^. Left: Mmin ~ 10^^ Mq 
and Right: deeper spectroscopic survey, M^i^ ~ 10^^ Mq. As the depth of spectroscopic survey increase, the constraints on b and G 
from the lensing-galaxy cross-correlation improves, because the stochasticity of the galaxy sample with respect to mass goes down when 
smaller halos are mapped. The constraint on fG also improves with the spectroscopic survey depth. The joint constraint on / and G 
also improves with spectroscopic depth at high z in particular. Adding lensing data to RSD splits the fG constraint into separate / and 
G constraints, which arc substantially more precise. 



the constraints become weaker at 2 > 0.5, indicating that 
the increasing shape noise and galaxy stochasticity dominate 
the improving mode counts. 

4.2 RSD constraint on fG 

Having RSD measurement alone, one can measure the pa- 
rameter fG after marginalizing over f/b. The green dashed 
lines in Figure [3] show the constraint on fG using the multi- 
tracer RSD method. We find fG is better constrained at 
high z, opposite to the lensing constraints on b and G shown 
in the previous subsection. The gain at high z for a{fG) 
mainly comes from having more modes as the survey volume 
dV/dz grows with redshift. When Mmin is smaller (compar- 
ing the right-hand panel to the left), a{fG) also drops, as 
expected, since we have more galaxies with a broader range 
of biases. 

BCll show, and we confirm here, that there is little 
change in the cosmological constraints from the multi-tracer 
RSD analysis from allowing the shot noise level to be a free 
parameter instead of fixing the Poisson value. 

4.3 Combined constraint on / and G 

The b and G measurement from lensing-galaxy cross- 
correlation in transverse modes can be added to the RSD 
analysis in the 3-D spectroscopic redshift survey over the 



same volume. This will help to break the f-b degeneracy 

existing in the case when RSD alone is available. Separate 

constraints on / and G are then achieved after marginaliza- 

tion over all the bias parameters, leaving 19 /'s and 19 P's 

~ f. p. 

in the Fisher matrix. An example of FLg|^g_,_jj^gj-, is shown 
in the right hand panel of Figure [T] As expected, the f — P 
Fisher matrix exhibits very little correlation between differ- 
ent redshifts. 

The solid lines of Figure |3] plot example constraints on 
/ and G vs 2 after marginalizing over all other parame- 
ters. The redshift dependence of (j(/) is similar to that of 
(t(/G), and the combined constraint on G improves signif- 
icantly over the case when lensing alone is available. This 
improvement is more prominent at high z. 



4.4 Constraint on gravity 

In the fiducial ACDM cosmology, a change in 7 results in 
predictable changes in the growth history (/ & G). There- 
fore, constraints on growth can be translated into a con- 
straint on 7. We present 17(7) as a function of Mens, the num- 
ber density of a lensing survey, and Mmin, which is equiva- 
lent to the survey depth of a spectroscopic redshift survey. In 
Figure |4l we compare four different cases: having RSD alone 
(I) (top left), lensing alone (II) (top right), RSD plus lensing 
over the same volume (III) (bottom left) and RSD plus lens- 
ing in separate volume (IV) (bottom right). All results we 
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Figure 4. Log of uncertainty (7(7) on growth parameter 7 as a function of spectroscopic survey limit A/min a-nd lensing source density 
Aliens- Color scale and contour levels are identical for all panels: Upper left: RSD measurement only; Upper right: lensing tomography 
only (including cross-correlation with photo-z galaxy samples); Lower left: lensing survey + galaxy redshift survey over the same volume 
at < 2 < 2, including cross-correlation between spectroscopic sample and lensing; Lower right: the same as lower-left panel, but the 
two surveys are not overlapping, so there is lensing cross-correlation with photo-z galaxies but not the spectroscopic sample. Note that 
in this case, the total number of modes are more than the case shown at lower left, i.e. the transverse modes of the two surveys are now 
independent. 
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Figure 5. Le/t; the square of the ratios of c(7) from RSD measurement versus that from lensing+RSD measurement. This is essentially 
showing how much larger a survey area the RSD measurement would need to achieve the same accuracy on 7 as when lensing information 
is added. Right: the same as the left but showing the ratio of lensing+RSD from separate volumes versus lensing+RSD over a common 
volume. The overlapping survey is most beneficial along a band with Af^in fa IO^^Mq. For deeper spectroscopic surveys, the separate 
surveys are nearly equivalent to overlapping, because the lensing information is a weak addition whether or not overlapping. Likewise for 
very shallow spectroscopic surveys and deep lensing survey (upper left corner), the lensing information dominates and overlap of RSD is 
irrelevant. 



show assume that the primordial CMB power spectrum is 
known, unless specified otherwise. Figure[5]plots a.y vs Mmin 
at two distinct Anions values, for different strengths of prior 
knowledge of the amplitude of Pcmb- Marginalization over 
the normalization of Pcmb leads to aj constraints about a 
factor 2 worse than knowing Pcmb exactly, but a prior with 
0.5% accuracy on Pcmb recovers almost all of this loss. 

4.4.1 (I) RSD alone 

When using multi-tracer RSD from a spectroscopic redshift 
survey, the result strongly depends on the survey depth, 
or Mmin (upper left plot of Figure [4] and dashed lines 
in Figures |6l By surveying halos with down to Afmin ~ 
RSD alone can already constrain 7 to ~ 
1%. The current GAM A survey reaches this survey depth 
iRobotham et aLll201ll ) up to z ~ 0.5, but one need to ex- 
pand the survey to cover half of the sky in order to achieve 
this accuracy. As Afmin increases, (7(7) increases rapidly — 
cr(7) ~ 2% for Mmin ~ lO"/i"^M0 and (7(7) ~ 10% for 
Mjnin ~ 10^^/i~^Mq (equivalent to galaxy number density 
of lO-'^-lO-'^Mpc"^). 

Central LRGs are considered as good samples for RSD 
measurement, as they reside at the center of their host halos 
theref ore should be free from the non-linear fin ger-of-God 
effect jOkumura fc Jing||201ll : lHikage et al.ll201ll '). They are 



hoste d by halos with Mmin > lO^^/i ^Mq (e.g. IZheng et al.1 
120091 ). LRGs are, however, an incomplete sampling of halos 
near 10^^ Mq. Therefore, RSD with LRG samples may 
not be as powerful as the Mmin = W^'^h~^ Mq forecast here. 



4.4-2 (II) lensing alone 

When a lensing survey alone is available, we can still esti- 
mate the mass density in the transverse modes using photo-z 
galaxies. In this case, the stochasticity between the galaxy 
density field and the mass field may be larger, so we conser- 
vatively assume End ~ 0.5 and add a weak prior on the A/" 
term as we have discussed in section [53] 

For this case of joint shear/density tomography without 
spectra, (7(7) varies by a factor of ~ 5 when Aliens varies 
from 1 to 100 (Figure 3] upper right). The upcoming DES 
is expected to have Mens ~ lOarcmin"^, which if scaled to 
/sky = 0.5 constrains 7 down to ~ 2% (Figure[6l left). Future 
surveys like LSST or Euclid, attaining A'lcns ~ 40 arcmin"^ 
in the most optimistic scenario, yield ~ 1.25 x reduction in 

Removing the CMB constraint on the amplitude of the 
primordial power spectrum may increase the error in 7 by a 
factor of 2. 
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Figure 6. Comparing <j{-y) as a function of Afniin from lonsing measurement (dotted lines), RSD measurement (dashed lines) and 
RSD+Lensing (solid lines) with different priors on PcMB in different colors. Left: the number density of galaxies in the lensing survey 
is Aligns = lOarcmin^^. Flight: ultra-deep lensing survey wtih Ni^^^^ = 100/, arcmin~^. When M^iin is large so the number density of 
halos in the galaxy redshift survey is small, constraints on 7 mainly come from lensing tomography. When the spectroscopic survey is 
deep, addition of lensing data does not improve the constraint on 7 substantially. If we had forecast a single-tracer RSD analysis instead 
of multi-tracer RSD, then the gain from adding lensing to a deep spectroscopic survey would be larger. Note that adding 0.5% prior 
on PcMB improves the constraint on 7 by a factor 2, and is close to the case of knowing Pcmb perfectly. The red stars indicate where 
RSD-|-Lensing is worse than lensing alone in the constraint of 7. This is because in the lensing alone case, we are using the projected 
galaxy-galaxy clustering from the lensing photo-z sample. The stochasticity in this case can be lower than that of the spec-z sample, 
when Muiin is very large. 



4-4-3 (III) RSD + lensing (same volume) 

Combining the two surveys will in general help to improve 
the constraint on 7. The cr(7) for the overlapping surveys 
is plotted in the lower left of Figure [H The left-hand plot 
of Figure [5] quantifies the improvement from adding lensing 
to the spectroscopic survey as the inverse square of the im- 
provement in (7(7), which is equivalent to asking what factor 
more survey area the RSD survey would require in order to 
match the improvement from the addition of lensing data. 
The amount of improvement depends on many factors. 

When the spectroscopic redshift survey is very deep, i.e. 
Mmin ~ 10^^/i"^Mo, RSD alone can already measure 7 at 
sub-percent level, if combined with 0.5% constraint on the 
primordial power spectrum from the CMB (Figure O. The 
number density of halo redshifts in such a survey is Hspec ~ 
10~^h^ Mpc""^, requiring a total of ~ 10* redshifts over half 
of the sky. In this regime, the improvement in the constraint 
of 7 by adding lensing data is very minor, and changing 
the lensing survey depth will not affect the result. This is 
consistent with the result of BCll (see their Figure 3). 

Figure [3] displays dramatic gains in constraint of / and 
Gat z > 1.2 when M^in is reduced from 10^^ to lO^^/i'^M©, 
yet only modest gains in (7(7) are seen in Figure [6] or on the 
left of Figure (5] This is because both the absolute values of 



dP/d'y and df /d'y become smaller at high 2, where is 
close to 1, so measures of / and P at z > 1 are less valuable 
in constraining gravity under this parameterization. 

When only halos with M > W^^h~^ Mq are targeted 
in the spectroscopy survey, the benefit of combining with 
a lensing survey becomes more prominent: equivalent to a 
factor of 2 to 3 increase in survey volume at Mmin ~ 3 x 
lO^^h'^MQ, and more than IQx when M^in ~ IO^/i'^Mq! 
(See left of Figure [S] and Figures [§1) However, even in this 
regime, lensing is never completely dominant in the range of 
Mens we are considering, in the sense that Figure |4] shows 
that improving the spectroscopic survey depth is always sub- 
stantially beneficial for measuring 7. 



4-4-4 (IV) RSD + lensing (separate volume) 

Here we compare the power of lensing and RSD surveys 
conducted over a shared /sky = 0.5 to surveys that do not 
overlap, covering distinct volumes. Having two surveys in 
separate volumes has the advantage of having twice as many 
transverse modes as the case of overlapping survey volumes; 
is this advantage outweighed by knowing the bias of the 
spectroscopic survey galaxies through the overlapped lensing 
survey? 
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To forecast the 7 constraints from separate surveys, we 
make the following alterations to the Fisher methodology 
for the combined surveys: first, we construct Flchs under 
the assumption that photo-z samples are being used for the 
galaxy density map (=> Ead = 0.5). Then we marginalize the 
bi values in FLens to leave constraints over only the Pi. For 
the RSD Fisher matrix, we marginalize over all the bia since 
no lensing constraints are available, leaving behind only a 
constraint on the product fP at each redshift. We also allow 
the RSD analysis to use all transverse modes, since these are 
no longer redundant with those in the lensing survey. The 
lensing and RSD Fisher matrices can again be summed and 
projected onto a single (7(7), plotted in the lower right of 
Figure |4] We also plot, on the right-hand side of Figure O 
the effective area gain of the overlapping survey relative to 
separate surveys. Note that this area "gain" could be as low 
as 0.5, i.e. a loss, since the combined survey does cover only 
half the volume of the separate ones. 

We find that having two surveys over the same volume 
is better than having them separated except for extremely 
deep lensing or spectroscopic surveys. The improvement is 
equivalent to a factor of 1.5 to 2 in survey volume when 
10^^ H'^Mq < Mmin < lO"/i"^M0, but very minor when 
the spectroscopic redshift survey gets deeper. This is true 
when the primordial CMB power spectrum is known to bet- 
ter than 0.5%. If we do not employ any CMB constraint 
and marginalize over Pcmb, then the 7 constraint will be 
degraded for each case, but the gain of having overlapping 
survey volume versus separate volume is larger, e.g. a factor 
of 3-4 in the regime when Mmin is large. The area gain fac- 
tor is ~ 1.5 even when the spec-z; is deep. Therefore, having 
a weak CMB prior makes the idea of combining two surveys 
over the same volume more useful, while a strong CMB prior 
help to reduce (7(7) in both cases and narrows the difference 
between them. 

Notice on the lower-right of Figure UJthat in the regime 
when Afmin > 10^^'^ Mq (shallow spectroscopic redshift 
survey), the constraint from lensing measurements is dom- 
inant and the depth of the non-overlapping spectroscopic 
survey becomes irrelevant. 

In summary, combining two surveys help most, rela- 
tive to separate surveys, when the spectroscopic redshift 
survey is modestly sparse, Mmin ~ W^'^'^ Mq, in the 
range of LRG surveys. When the spectroscopic survey is 
deep (Mmin < W^^h~^ Mq), it dominates the error budget 
and it matters less whether the lensing survey overlaps or 
not. On the other hand, when the spectroscopic survey is 
very shallow (A/min > lO^"*-^/i~^M0), then even a modest 
lensing survey (A'lens > 5arcmin~'^) dominates the infor- 
mation, and it matters less whether the spectroscopic sur- 
vey is coincident. There is, however, no regime of feasible 
large-scale surveys for which the separate surveys constrain 
7 better than overlapping surveys. 



5 CONCLUSION AND DISCUSSION 

We have shown from Fisher matrix forecast that the con- 
straint on the growth of structure and gravity can be reduced 
percent-level or even sub-percent level by combining a spec- 
troscopic redshift survey with a photo-z weak lensing survey 
over the same volume. Whereas BCll merely assumed that 



some measure of galaxy bias was available to add to RSD in- 
formation, we verify here that a realistic tomographic weak 
lensing survey does in fact yield bias information sufficient 
to realize a substantial gain in accuracy on the growth pa- 
rameter 7. 

Following the sugg estions of [Pe3 (|2004l ') we use the 
shear-galaxy cross-correlation to measure the galaxy bias in 
the transverse modes — a measurement which is free of sam- 
ple variance — and apply it to the multi-tracer RSD analysis 
in a spectroscopic reds hift survey (|McDonald fc Seliakll2009l : 
iBernstein fc Caill201ll '). The combination of the two surveys 
make it possible to measure the linear growth function G 
separately from its derivative / = c?lnG/dlna, whereas 
RSD alone can only measure the product fG. 

The performance of multi-tracer RSD measurement de- 
pends on the spectroscopic survey depth, the range of galaxy 
biases in the sample, and the number of linear modes avail- 
able. The performance of the shear -|-galaxy analysis on the 
transverse modes depends on: (1) the level of stochasticity 
between the galaxies and the projected mass, (2) the depth 
of the lensing survey, or shape noise. When combining two 
measurements together over the same volume, the results 
will depend on all those factors that affect each of the sur- 
vey. 

We have demonstrated that for the constraint of the 
7 parameter, combining two surveys is better than having 
each of them alone, roughly a factor 1.5 improvement (in 
survey-area terms) in the regime of likely feasible surveys: 
source density Mens ~ lOarcmin"^ in the lensing survey, 
and galaxy surveys complete for halos in the cluster or small- 
group range Mmin = 10^'^-10^*/i~^A/q, similar to LRG sur- 
veys. For A/niin > lO^^/i~^Af0, the lensing-fRSD survey has 
constraints many times more powerful than the RSD survey 
alone. The 7 parameterization of growth predicts very little 
change at 2 > 1; a different model for deviations from Gen- 
eral Relativity could gain even more from the combination 
of lensing and RSD surveys. 

Having prior constraints on the amplitude of the pri- 
mordial power spectrum from the CMB is useful in general. 
Knowing Pcmb to 0.5%, easily within the statistical power 
of Planck, garners most of the « 2x gain in accuracy on 7 
that is possible with perfect a priori knowledge of Pcmb . If 
Pcmb is more poorly known, the gain of having overlapping 
surveys over the case of separate survey vo lume is increased. 

During preparation of this paper, iGaztanaga et all 

(|201lD released very similar calculations of the benefit of co- 
incident lensing and spectroscopic surveys. Their assumed 
survey configurations and free parameterizations differ sub- 
stantially from ours, so direct quantitative comparison is 
not possible. In the particular case of constraints on 7, they 
find overlapping surveys reducing (7(7) by ~ 2.4 x compared 
to seperated survey volumes, equivalent to a 6x increase 
in survey area in the language of our Figure [5] where we 
find ~ 1.5 X areal gain. This is qualitatively consistent with 
our conclusion, but the origin of the substantial quantita- 
tive difference is difficult to ascertain gi ven the different 
assumptions about survey characteristics. IGaztanaga et all 
(2011.) also find substantial gains in accuracy of dark energy 
equation-of-state determination from overlapping surveys. 
Our analysis holds this fixed so we would not have detected 
these gains; we plan to broaden our analysis to the case of 
unknown distance-redshift relations in the near future. 
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We notice that bias measurement can in principle also 
be measured using the same spectroscop ic sample from the 
galaxy bispectrum (jSimpson et al.l^2011^ . If the same accu- 
racy of bias can be obtained in this method as using lensing, 
one can simply use one spectroscopic redshift survey to ob- 
tain the same measurement, which might be another attrac- 
tive survey strategy since no lensing survey is needed. The 
lensing survey is, however, a straightforward measure of the 
galaxy bias, free of assumptions about perturbation theory, 
second-order bias, and other issues with the bispectrum. 

Use of smaller-scale modes are attractive in the sense 
that one may gain many more modes from the same volume 
of survey. Growth test statistical accuracy improves rapidly 
with increasing fcmax- Non- linear effects in the density or 
velocity field and scale-dependent bias may, however, ruin 
the attempt to achieve percent-level constraint on param- 
eters. Efforts have been made to improve RSD predictions 
for smaller-scale modes |Scoccimarrol 20041 : Ijennings et al.l 
I2OIII : iHikage et al.ll201lljTang et al.ll201lh . e.g. fcmax ~ 0.3, 
though it is important that predictions be made for galaxies 
or ha los rather than all m ass parti cles in an A'^-body sim ula- 
tion ijjennings et al.ll201ir i. e.g see lReid fc White! (|201ll l for 
modeling of halos. Better understanding of the non-linear 
biases of different tracers is required before one can confi- 
dently select the fcmax that admits the most accurate growth 
constraints. 
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APPENDIX A: OPTIMAL WEIGHTING AND 
MINIMAL STOCHASTICITY 

We take the halo model description of the stochasticity E 
that has been developed in CBS. The basic idea is to split 
halos into different mass bins, apply the optimal weight uiopt 
to each of them: 
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where 
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where the first term on the right hand is equivalent to mass 
weighting, and the second term being close to bias weighting, 
^(m) is the halo bias respect to the 'continuous halo field', 
u{k\m) is the Fourier transform of the NFW halo profile 
(NFW), p is the mean mass density. 
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We then obtain the corresponding stochasticity between the 
weighted halos field and the mass field: 
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This model has been shown to be in good agreement with 
simulations (|Cai et al.ll201ll ). We will use -E^opt BuS the fidu- 
cial value of stochasticity in this paper. 



